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of the reinforcement from 94% to 162% of the volume of the 
material removed by the hole, the number of reinforcing 
layers (1 to 4) and the outer radius of the reinforcement 
layers from 0.75 inches to 1.5 inches. The orientation of 
the reinforcing layers was either 45° or 0° to the load 
direction. Results of the computational analysis indicated 
that the reduction of maximum strain (in the direction of 
the load at a point at the edge of the hole 90° from the 
load direction) by the reinforcements was at most 12%, with 
apparent dependence on the number of layers used and the 
volume of the reinforcement. Experimental testing of three 
of the configurations confirmed the accuracy of the finite 
element analysis and demonstrated that the reinforced panels 
recovered 5 to 10% of the basic laminate strength above that 
of the reinforced panel. In addition to the computational 
analysis, two isotropic empirical predictions for stress 
(strain) concentration at the edge of a circular hole were 
examined to determine their applicability to asymmetric 
reinforcement of a quasi-isotropic panel, with mixed results. 
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ABSTRACT 


The response to uniaxial loading of a graphite epoxy 
panel containing a one-inch diameter hole was analyzed by 
finite element analysis and tested experimentally. The 
analysis modeled the basic unreinforced panel and six dif- 
ferent asymmetric reinforcements consisting of additional 
layers of the same material of circular shape on one side 
of the panel laminate. The reinforcement configurations 
varied the volume of the reinforcement from 94% to 162% of 
the volume of the material removed by the hole, the number 
of reinforcing layers (l to 4) and the outer radius of the 
reinforcement layers from 0.75 inches to 1.5 inches. The 
orientation of the reinforcing layers was either 45° or O= 
to the load direction. Results of the computational analysis 
indicated that the reduction of maximum strain (in the direc- 
tion of the load at a point at the edge of the hole 90° from 
the load direction) by the reinforcements was at most 122%, 
with apparent dependence on the number of layers used and 
the volume of the reinforcement. Experimental testing of 
three of the configurations confirmed the accuracy of the 
finite element analysis and demonstrated that the reinforced 
panels recovered 5 to 10% of the basic laminate strength 
above that of the reinforced panel. In addition to the 
computational analysis, two isotropic empirical predictions 
for stress (strain) concentration at the edge of a circular 
hole were examined to determine their applicability to 
asymmetric reinforcement of a quaSi-isotropic panel, with 


mixed results. 





tT 


if boil 


IV. 


VI 


TABLE OF CONTENTS 


INTRODUCTION ---- 32-29 2-2 = 
COMPUTATIONAL ANALYSIS --c3 rococo oo enn 
A. MODEL FORMULATION ----------------------- 
B, COMPUTATIONAL RESULTS ---r-c-o-n orn 
1. Computed Surface Strains ------------ 
2. Contour Plots ----------------------- 


3. Computed Strain Concentration 


Factors -------<<-------------- ------ 
ALTERNATE THEORETICAL ANALYSIS ---2---------- 
A. UNREINFORCED PANEL ----<2-23-3-2e3%<2--------- 
B. REINFORCED PANELS -----"-----2------------- 

i. “Empirical Isotropic Approaches -=--—-= 
EXPERIMENTAL METHOD ---9--- 9-9-9" ---------- 
Pe eeannb SHEEECTION AND CONSTRUCTION -=------- 
B. TEST APPARATUS rrr rrrettte treet err sr----n 
C. TEST PROCEDURE ---9-2---- 02-22-2205 ------ 
EXPERIMENTAL RESULTS -----------2"--20-------- 
A. UNREINFORCED PANEL ---------------------- 


B. REINFORCED PANELS wrrrrtrrrterrrtrr cre 


{PE sPERIMENTAL DETERMINATION OF K,, “===<---- 


90 
DISCUSSION OF RESULTS -<-<-----<-<-<<---------- 
A. COMPUTATIONAL RESULTS ------------------- 
B. EXPERIMENTAL RESULTS -------------------- 
C. COMPARISON OF PREDICTED AND MEASURED 
VALUES OF Rog = eS ee eee 


16 
At 
Zi 
29 
29 


36 


40 
48 
48 
49 
49 
25) 
oe 
58 
62 
64 
64 
65 
dal. 
78 
78 


82 


88 





WMElE 0) ee 2 ee = 23 


VIII RECOMMENDATIONS --------------------------------- 95 
APPENDIX A NONREINFORCED PANEL----------------------- 96 
APPENDIX B FIRST REINFORCEMENT CONFIGURATION -------- 108 
APPENDIX C SECOND REINFORCEMENT CONPIGURATION ------- 126 
APPENDIX D THIRD REINFORCEMENT CONFIGURATION -------- 136 
APPENDIX E FOURTH REINFORCEMENT CONFIGURATION ------- 156 
APPENDIX F FIFTH REINFORCEMENT CONFIGURATION -------- 17 
APPENDIX G SIXTH REINFORCEMENT CONFIGURATION -------- 184 
LIST OF REFERENCES ------------------------------------ 205 
INITIAL DISTRIBUTION LIST ------------------------~---- 206 





jgi8 


TL 


IV 


Wa 


el 1 


at T 


IX 


XI 


LIST OF TABLES 


COMPUTATIONAL MODEL CONFIGURATION ------- 


COMPUTED STRAIN CONCENTRATION 
FACTORS wre rrr nn nn 


Poo TROLrTe. - REDErCTION oF K 99 se eam ee 
ALTERNATE ISOTROPIC PREDICTION OF Kgo7777 
MATERIAL PROPERTIES OF HMF 330/34 ------- 
QUASI-ISOTROPIC LAMINATE PROPERTIES ----- 
EXPERIMENTAL FINDINGS ------------------- 
EXPERIMENTALLY COMPUTED VALUES OF KgQ7-77 


STRAIN CONCENTRATION FACTOR BASED 
ON MAXIMUM STRAIN ----------------------- 


COMPARISON OF EXPERIMENTAL AND 
COMPUTED MICROSTRAINS 3------------------- 


COMPARISON OF MEASURED AND PREDICTED 
VALUES OF Kgo orn nnn = = 


27 


43 


oy 


De 


24 





10 


Jug 


2 


3 


14 


aL} 


16 


iL / 


18 


IY, 


20 


iol were LGORES 


Mesh Generated to Model Panel ------------------ 
Cross Section of a Typical Reinforcement ------- 
Comparison of Front Surface Strains (€_) ------- 


Comparison of Front Surface Strains (e,) ------- 


Koq 


K, Vs Number of Reinforcement Layers ----------- 


0 


Koq 


Test Panel Dimensions ----<---------------------- 
Experimental Apparatus ------------------------- 


Strain Gauge Placement on All Panels — 


Comparison of Computed and Measured Strains 


EC. Vs X/R Basic Panel -------------------------- 
Comparison of Computed and Measured Strains 

€e Vs X/R Basic Panel -------------------------- 
Comparison of Computed and Measured Strains 

Ee, Vs Z/R Basic Panel ------------------+-------+- 
Comparison of Front and Back Strain 90° From 
Load Axis - Basic Panel Gauges 5 and 19 -------- 
Typical Panel Fracture -------------<+----------- 


Measured Tangential Strains at the Hole Edge 


for a Load of 30,000 Lbs ---------------------+-- 
Strain Ee, Vs Load at Hole Edge ---------------- 
Nonlinear Behavior of Panel 1B 

EC, Vs X/R at the Hole Edge --------------------- 
Comparison of Isotropic Prediction and 

Measured Values of Kgq NN eon ae eas ol ie = 


Comparison of Isotropic Prediction and 


Measured Values of Kgq oe ee oy ee ee 


Vs Number of Reinforcement Layers ---------- 


Vs Reinforcement Volume -<-<--<-<-----—----------- 


ZS 


Sy) 


38 


45 


46 


47 


56 


ao 


61 


66 


67 


68 


69 


Vek 


ig 


74 


87 


jek 


92 





Ze 
22 
23 
24 
25 
26 
27 
28 
29 
30 


31 


52 


33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 


45 


Computed ¢«_ Vs X/R Unreinforced Panel ------------ 96 


Zz 
Computed ¢, Vs X/R Unreinforced Panel ------------ 97 
Computed e,, Vs X/R Unretntorced Panel =--<--------- 98 
Computed e, Vs 2/R Unreinforced Panel ------------ 99 
Computed ¢,, Vs Z/R Unreinforced Pane = ———— 100 
Computed e, Vs 2/R Unreinforced Panel ------------ abel 
Contour Plot IS a a LO 2 
Sontour Plot of ee IMA ae= SS oo 103 
SOnmtour Or No ee ee ee 104 
Contour of Out of Plane Deflections -------------- 105 
Comparison of Measured Front and Back Strains at 
Saige oel> amamc0e— basic Panel) =—=—=—=<-—--—---—-----— 106 
Comparison of Measured Front and Back Strains at 
Gauges 9 and 21 —- Basic Panel -----<-<-<<<<<-<------ 107 
aS Vs X/R erm wer rr tr rrr rrr et rest sressce-- 108 
on Vs X/R seer rt rn rr rer 109 
— VS X/R srrer ttre sn tr rr rr rr otter cece IIL 
E. VS 2/R erm rrr terre rr rrr rr rr ere ital 
€ Vs 2/R were rest tt rr rr renee eZ 
— Vs Z/R wre etr wrt sett rr rrr rrr NES: 
Contour Plot of No Configuration lA -------------- 114 
Contour Plot of Ny Configuration 1A -------------- pas 
Eontour PLOt of Nye Configuration 1A ------------- 116 
Contour Plot of N. Configuration 1B -----22r-e+-+----- Ibley 
Contour Plot of N. Configuration 1B -------------- 118 
Contour of Nos Configuration 1B ------------------ 119 


Comparison of Measured Front and Back Strains at 
Gauges 5 and 19 ~- Configuration 1B ----<s-<-------- 120 





46 


47 


48 


49 


50 


al 
2 
De 
54 
=) 
56 
5 / 
58 
59 


60 


61 
a2 
63 
64 
65 
66 
67 


68 


Comparison of Measured Front and Back Strains 
Macc omeoeandeay — COntighuration 1B ----------- 
Comparison of Measured Front and Back Strains 
at Gauges 9 and 21 =- Configuration 1B -e--------- 


Comparison of Computed and Measured Ee. Vs X/R 
Configuration 1B --------------------=------------ 


Comparison of Computed and Measured Strain 
ee OR Oot land Benge pam ao a a n= 


Comparison of Measured and Computed Strains 
Bop WS AY MR COMES USAGES seal NN aa 


Gontour Plot of Ny Configuration 2A -------------- 


Contour Plot of N,, Configuration 2A ------------- 


z 


Contour of Out of Plane Deflections 
Configuration 2A ------------------------------- == 


1d 





69 


70 


yal 
2 
73 


74 


a5 


76 


ia 


78 


no 


80 


81 
82 
83 
84 
85 
86 
S7 
88 
89 


90 


Contour of N,, Configuration 3A ----------------- 


Contour of Out of Plane Deflections 
Configuration 3A --------- 9-2 ee 


Conmeour Of N, Configuration 3B ------------------ 
Contour of Ny Configuration 3B ------------------ 
Contcurellor o£ Nye Configuration 3B ------------ 


Contour Plot of Out of Plane Deflections - 
Configuration 3B -------------------------------- 


Comparison of Measured and Computed Strains 


a Vs X/R - Configuration 3B ------------ 9-H 
Comparison of Measured and Computed Strains 

Ey Vs X/R Configuration 3B -------- 
Comparison of Measured and Computed Strain 

ec, Vs Z/R - Configuration 3B -------------------- 
Comparison of Measured Front and Back Strains 

at Gauges 5 and 19 - Configuration 3B -----<------ 
Comparison of Measured Front and Back Strains 

at Gauges 15 and 19 - Configuration 3B ---------- 
Nonlinear Behavior of Panel 3B 

Ee, Vs X/R at the Hole Edge ----------------=------ 
c. VS X/R wren rr = 
ce, Vs X/R wee 
ce VS X/R wer ee = 
E Vs 2/R were + 
E Vs 2/R eee re == 
2a VS 2/R eer 
Gomtour PLOt of N, Configuration 4A ------------- 


Contour Plot of Ny - Configuration 4A ----------- 
Contour PLOt Of Nye - Configuration 4A ---------- 


Contour Plot of Out of Plane Deflections 
Configuration 4A -------------------------------- 





gL 
a2 
23 


94 


95 
96 
oi 
98 
29 
100 
101 
102 
FO 3 


104 


105 
106 
uO7 


108 


109 
LO 
Jal 
ie 2 
£13 
114 
115 


COMEOuUL@E LOE Of N, - Configuration 4B ------------ 166 


Comeour PLot or Ny - Configuration 4B ------------ iG? 
Contour of See - Configuration 4B ---------------- 168 
Contour Plot of Out of Plane Deflections - 
Configuration 4B qe e nn ne ene eres sere--- 169 
€ Vs X/R -oo-- 2-2 -  -  -  -- ---- ----- 170 
— VS X/R errr rrr rrr rr rrr rr rrr rennin e--- 171 
See VS XJR qe rr rr rr rrr rrr rere i 2 
B, WS BPR SESS la eae ata 7s 
E Vs Z2/R weer eee re 174 
ee = ----=-—— 175 
Geneour Plot of Nw - Configuration 5A ------------ 76 
Conecun Plot Of Ny - Configuration 5A ------------ 1g 
Contour Plot of ee Configuration 5A ----------- 178 
Soueolm Plot Of Out of Plane Deflections 
CS UESE TGA GENE TS RON TUE a a a hoe 179 
Sourotee blo t@iOine Na CONLIgGuraAELON 5B ----—----=-— 180 
GConteour Plot of ro Configuration 5B ------------ 181 
Concour Plot Of es - Configuration 5B ----------- iow 
Contour Plot of Out of Plane Deflections - 
Configuration 5B ---- 2-22-28 ne een reese 133 
ee ee 184 
SWS BYR Sear cea 185 
SRS = C7 aT) ec a a a ca a es se ee ee 

oe Vs X/R 186 
a Vs 2/R errr rrr tr rrr rr rn eee ey 
oS VS 2/R rrr rr rrr rr rrr rr rrr 188 
€ VS 2/ Ree rn 2 a a  - 189 
XZ 

eontourn Plot of Ne Configuration 6A ------------ 130 


eZ 





6 
ie 7 


118 


9 
120 
2 1 


eZ 2 


2 3 


124 


25 


#26 


2 7 


2 8 


29 


Concour Plot o£ Ny - Configuration 6A ---------- deal 


Contour Plot of Noe - Configuration 6A --------- 192 
Contour Plot of Out of Plane Deflections 

Configuration 6A weer rrr eet nr rr rr errr nn fern Eos 
Contour Plot of No - Configuration 6B ---------- 194 
Contour Plot of Neo Configuration 6B ---------- 195 
cancout BGt Of Noe - Configuration 6B --------- 196 


Contour Plot of Out of Plane Deflections - 
Configuration 6B nent ere ee eer rn rrr rrr nr 197 


Comparison of Measured and Computed Strain - 
Ee, Vs X/R Configuration 6B -----------------=---- 198 


Comparison of Measured and Computed Strain - 
Ee, Vs X/R - Configuration 6B ---~---------------- 199 


Comparison of Measured and Computed Strain 
Ey Vs Z/R Configuration 6B ----------- - 200 


Comparison of Measured Front and Back 
Strain at Gauges 5 and 19 - Configuration 6B-- 201 


Comparison of Measured Front and Back Strains 
at Gauges 15 and 20 - Configuration 6B --------- 202 


Comparison of Measured and Computed Strains 
at Gauges 9 and 21 - Configuration 6B ---------- 203 


Nonlinear Behavior of Panel 6B 
ce, Vs X/R at Hole Edge - Configuration 6B ------ 204 


IS 





CJ 


a = a 2 


v9) 


ba 2, XZ 


90 


LIST OF SYMBOLS 


cross sectional area of the reinforcement 
semi-major axis length of an elliptical hole 
semi-minor axis length of an elliptical hole 
Young's Modulus 

thickness of both reinforcement and panel 
Curvature 

stress moment 

stress resultant 

radius of hole, Q.5 inches 

thickness of panel 

volume 


coordinate axis, horizontal through plane of 
panel 


coordinate axis, through thickness of panel 


coordinate axis, vertical through plane of 
panel 


strain 

microstrain 

midplane 

components in plane of panel 


O 


O- to Load axis 


90° to load axis 


14 





ACKNOWLEDGEMENTS 


The author wishes to express his most grateful 
appreciation to Dr. Milton Bank for his enthusiastic 
Seadance, Dr. J. A. Bailie for his constant untiring 
support, Gordon Ferguson, Michael Murphy and Norman Cyr 
for their patient and detailed guidance through the DIAL 
COMPUTER CODE. In addition the author extends heartfelt 
thanks to the individuals at Lockheed Missiles and Space 
Co. for producing the test articles. Last, but certainly 
not least, the author offers many thanks to Mr. Glen 
Middleton and Mr. Bob Besel for their exhaustive assistance 


in setting up the test apparatus. 


5 





I. INTRODUCTION 


The steady increase in the use of composite materials in 
aircraft design has been made possible by years of basic 
research concerning the qualities of composites. Examination 
of the behavior of composites in basic design configurations 
has provided the answers required to use composites in those 
configurations with a high degree of confidence, leading to 
the incorporation of those configurations in evolving designs. 
Present applications of composites in design are quite natu- 
rally dominated by configurations that have been understood 
for the longest time, and have avoided where possible the use 
of configurations which have not been as deeply researched. 
Because the researchers in composite materials must pay care- 
ful attention to the effects resulting from the inherent 
orthotropicity of the material, such as stacking sequence of 
layers, laminate direction and composition, the time required 
to understand the behavior of composites as fully as we 
understand that of isotropic materials such as steel alloys 
is much greater. There are many areas in which basic research 
in composites is still in progress. 

One such area of research and possible application is 
that of a composite panel under tensile loading which contains 
an intentionally placed stress concentration such as a circular 
hole. A wing panel with an access hole or a tank having a 


hole in its wall to permit the attachment of a feed line 


aG 





would be examples of such an application. The interlaminar 
shear present near the hole edge, coupled with the stress 
concentrations created by the existence of the hole, pose 
complications which seem to be avoided whenever possible by 
design engineers, causing the components to be potentially 
lacking in flexibility of design. 

The response of composite panels containing stress concen- 
trations to uniaxial tensile loading has been the subject of 
past research. Rowlands, Daniel and Whiteside [Reference 1] 
determined that the material properties of the laminate 
specify the strain behavior in the vicinity of the hole (see 
Analysis, page 49). Also, there was evidence that inter- 
laminar stress did influence the response of the panel. The 
same work also demonstrated that a well constructed finite 
element model could be used to compute the behavior of the 
plate with a high degree of agreement with experimental 
results. 

The same research team [Ref. 2] also pursued the influence 
of laminate orientation and stacking sequence on the stress 
and strain behavior of a plate containing a hole. Their 
findings show that the orientation of the plies in relation 
Meche load axis (specifically the number of + 45° plies) 
has a significant effect on the reduction of tangential 
strain around the edge of the hole. They also showed that 
the reduction of the total strength of the panel, compared 
With an identical panel without a hole, can be mitigated by 


the predominant use of + 45° plies. In addition, the mode 


by 





of failure of the panel was shown to be dependent upon the 
stacking orientation and sequence. 

While the reinforcement of holes in composite panels 
would seem to be a logical extension of this line of research, 
little information concerning reinforcements exists in the 
published literature. McKenzie [Ref. 3] demonstrated some 
success with hoop wound graphite epoxy rings attached asym- 
metrically to both aluminum and graphite epoxy panels. The 
effect of varying the outer radius and thickness of the 
reinforcements was investigated, and the potential difficulty 
of reinforcing only one side of the panel was mentioned. The 
method used in the investigation had the disadvantages of 
separate production of the reinforcement and panel, attach- 
ment processing, and post attachment adjustment if dimensional 
consistency was required. 

The investigation reported here. studied the effectiveness 
of asymmetric reinforcement of a quaSi-isotropic graphite- 
epoxy plate containing a circular hole. Asymmetric reinforce- 
ments were chosen because of their greater potential ease of 
manufacture and their applicability in some instances, 1.ée., 
where mating of two panels would be adversely affected by 
small raised areas of reinforcement. The reinforcement 
Meeposea for study consisting of additional laminate added 
asymmetrically, i.e. to only one surface of the panel, and 
covering only a portion of the panel in the region of the 


hole. The reinforcements were arbitrarily chosen to be of 
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Circular shape arranged symmetrically about the hole and 
varying in radius to gradually butld up thickness near the 
hole. This choice was based on the historical use of such 
reinforcements in isotropic materials. In being laid up 
with additional layers of the same material, the reinforce- 
ments have the advantage of co-curing with the panel, with 
no thermal property variations and no complication on the 
underside tooling (which would be necessary if the reinforce- 
ment were symmetric). 

This investigation selected six different reinforcement 
configurations based on volume of the reinforcement, number 
of layers of reinforcement and outer ring radius. These six 
configurations were then modelled for finite element analysis 
to compute stresses, strains and displacements. In addition, 
because of the quaSi-isotropic nature of the panels, the 
parameters of the reinforcement configurations were used in 
two isotropic empirical analyses for comparison to the finite 
element results. Experiments were conducted on the quasi- 
isotropic nonreinforced panels with and without holes (to 
establish a baseline for comparison) and on panels having 
several of the reinforcement configurations. The experimental 
results were then compared with the finite element results to 
verify the models and to the alternate analysis to determine 
the accuracy of the more simplified approaches. The attention 
of the computations and comparisons concentrated on the maximum 


Strain of the panel, located at the edge of the hole 90° from 


JES, 





the direction of the load, based on previous studies of 
composite panels which indicated that the strain in this 
area was the primary cause of failure. Finally, the 
experimental and computational results were compared in 
order to draw final conclusions about the effectiveness 


of the reinforcements studied. 
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IIT. COMPUTATIONAL ANALYSIS 


A. MODEL FORMULATION 

The analysis of the basic panel and all reinforcement 
configurations was completed through the use of a finite 
element computer code known as DIAL [Ref. 4], developed by 
Lockheed Missiles and Space Co. This code is a series of 
programs which compile files in a data base, develop the 
mesh, loading and boundary conditions of the model, and 
finally compute the nodal deflections and forces. All 
files created in the running of the several interrelated 
programs are then used in the post-processing portion of 
the code, which creates the output data in terms of strain 
and stress resultants. The actual solution to the matrix 
equation is computed using a modified Newton-Rapson Eee manne 
and strictly linear analysis. 

Since the computational analysis was to predict the 
behavior of panels which was to be experimentally verified, 
the physical dimensions of the computer model were determined 
by the design requirements for the test panels. For the 
discussion concerning the choice of the panel size, see the 
section on Panel Selection and Construction. The model 
created to duplicate the actual panels used in the experiment 
took advantage of the symmetry of the problem to reduce the 
number of elements and degrees of freedom. One quarter of 
the plate was used, corresponding to the upper right hand 


Ze 





corner of the actual plate as normally viewed. The mesh. 

was generated by a standard mapping from the k-j plane to the 
x~z (actual) plane. Because of the strains and stresses were 
expected to vary more rapidly near the hole edge and in 
anticipation of the circular nature of the proposed reinforce- 
ment configurations, several closely spaced rings were placed 
meer che hole at radii of 0.625, Q.75, 0.875, 1.00, and 1.25 
inches. These dimensions were chosen to produce a sufficient 
number of elements near the hole for accuracy while keeping 
the total number of elements low enough to prevent excessive 
computational run time. The mesh generated by the computer 
is shown in Figure l. For the basic panel, the mesh uses a 
total of 262 elements and 853 node points, of which 45 
elements are contained within the 1.25 inch radius. 

The elements selected to occupy the mesh for the panel 
were standard isoparametric quadrilateral elements having 
e€ight side nodes and one interior node. The elements utilize 
two Gauss quadrature integration points per side and include 
two points through the direction of thickness, giving them 
the properties of a modified thick shell element. The mate- 
rial property definition of the code permitted the layering 
of individual laminae within the element by defining the 
material properties and orientation of the laminae, so that 
in effect, each element of the basic panel has the identical 
properties of the laminate it is to model. The code also 


made provision for the buildup of laminae over a portion of 
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Pigure l. 





\ Nee 


Mesh Generated To Model Panel 
(Perspective View). 
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the panel to represent the areas of reinforcement, with the 
same ability to define the material properties of the partial 
laminae. This in effect adds additional elements to the panel 
Eeyeplacing additional elements over the basic panel, but does 
not add any degrees of freedom. This provision restricts the 
placement of the reinforcement layers to the radii previously 
described above, but allows for one. reinforcement layer to 

lay upon a layer below it. An example of the way the panel 
and reinforcements were built up by this process is shown in 
mieure 2. 

The quarter panel mesh was constrained by symmetry boundary 
conditions along the lower edge and the right hand vertical 
edge as seen in Figure l, two lines of symmetry passing through 
the center of the circular hole. In addition, to model 
tensile loading the upper edge was permitted to move in the 
plane (x and z deflections) and no rotations were permitted 
along that edge, which corresponds to the portion of the panel 
just at the edge of the load bearing fiberglass tabbing. All 
other nodal degrees of freedom were automatically made con- 
Sistent by the code. By the application of these boundary 
conditions, the problem was reduced to one having a total of 
4824 degrees of freedom in the solution. 

The loading of the panel was modeled as a uniformly 
distributed load of 1000 pounds per inch across the top edge 
of the panel, in the positive Z direction. This corresponds 
tO a total tensile loading of 10,000 pounds. The distributed 


loading was appropriately factored into concentrated loads at 
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Figure 2. Cross Section of a Typical Reinforcement. 


Zz 





each of the 29 nodes along the top edge. Because of the 
linear solution technique, results of this one loading can 
be adjusted to show the effects of any similar concentrated 
loading by simply multiplying the results by the ratio of 
the other load to 10,000 pounds. 

As previously discussed, the reinforcements chosen for 
Study were circular in nature and approximately equal in 
volume to the material removed from the panel by the hole. 
The six reinforcement configurations discussed were also 
divided into two subconfigurations, representing: a) all 
faminae oriented + 45°; or, b) one or two laminae closest 
to the panel having an orientation of q°. Table I shows 
the different reinforcement configurations and the labels by 
which they will hereafter be referred. 

The choice of 0.75 inches for the minimum outer ring 
radius was based on the concern that smaller reinforcements 
would be less effective due to the relative size of the 
region near the free edge affected by interlaminar stresses. 
Since the interlaminar stresses create a three dimensional 
Stress state in a region about one thickness from the hole 
edge (a total radius of approximately 0.6 inches for these 
panels) [Ref. 5], the minimum outer radius should exceed that 
distance by a reasonable amount in order to be effective. 
(The variation of ring radii in Q.25 inch increments was 
partially dictated by the mesh generation and also because 
it created an effective ramp between the layers of approxi- 


mately 3° for Single layered rings and 6° for double layered 
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TABEE 


COMPUTATIONAL MODEL CONFIGURATIONS 


CONFIGURATION RING THICKNESS ORIENTATION 
2 VOL RADIUS, IN DEG 


LA ‘ 3026 +45 
162% ee O02 6S -45 
1B 0.75 0.028 +45 
162% 1:25 O.0:2.8 0 
2A 1.00 0.028 +45 
103% d bap a 0.028 -45 
2B ire0'0 Q.028 +45 
103% ie 5 O10) ars: 0 
3A O75 Q.014 +45 
1183 1.00 0.014 -45 
eee 0.014 +45 
3B Oi75 0.014 +45 
1183 ie O0 0.014 -45 
IEA a 0.014 0 
4A G29 5 0.028 +45 
96.9% 1s eS) 0.014 -45 
4B Oras 0.028 +45 
96.9% Te25 Oeul4 Q 
5A 0.75 0.028 +45 
106.2 SG 20125 -45 
5B 2275 0.028 +45 
106.2 OO 0.028 0 
6A 15 0 0.014 +45 
100.0 
6B 155 0 Q.014 0 
100.0 
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rings.) These small ramp angles from one ring to the next 
were intended to provide an easy load path for the stresses 
in the panel, thus reducing stress concentrations in the : 
transition region from one ring to the next. 
The choices of the number of layers and radii were direc- 
ted by two considerations: a) keeping the total volume of 
the reinforcements close to that of the material removed by 
the hole, to reduce the weight added to the structure and 
also to permit comparison to other investigations; and, b) 
avoiding configurations which are unusually thick in the 
vicinity of the hole, resulting in extremely high bending 
loads and potential problems in mating with adjacent surfaces. 
For most of the twelve subconfigurations and the basic 
panel without reinforcement, the following results were 
obtained: 


1. Contour Plots of the Stress Resultants 


a. No 
lo MY 
x 
GaN 
xz 


dad. Out of plane deflections (not available for 
Hilscontigurat tons) 
2. Element strains and curvature for elements on the 
two lines of symmetry. 
3. Element Stress Resultants and Stress Moments (N's 


and M's) for the same elements as above. 
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For convenience, all of the figures for the basic panel 
and for each reinforcement configuration are separated into 
individual sections found in Appendices A-G. The results 


contained in the Appendices are discussed below. 


B. COMPUTATIONAL RESULTS 
1. Computed Surface Strains 

Because the experimental portion of this investiga- 
tion used strain measurements to determine the response of 
the tested panels, the numerical output from the DIAL program 
for all of the configurations was used to compute the value 
of surface strains on the front (reinforced side) and back 
of the panel. The values of the panel midplane strains and 
curvatures were used to compute the surface strains using the 


classical laminate theory strain equation 


€ € k 
4 x >-& 
e — 
exz xz ee 
: J 5 , 


The values of the surface strains were computed by employing 
the output of midplane strains and curvatures of one integra- 
tion point per element to reduce the amount of data reduction. 
The integration point selected was the point closest to the 
edge of the hole and axis of symmetry. For consistency, the 
same integration point was used for all other elements. In 
this way, a total of 11 points was used for calculations 


along the X axis (across the direction of the tensile stress) 
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and 16 along the Z axis (parallel to the stress). The 
results of these calculations are found in the first six 
figures of each Appendix, with both subconfigurations shown 
on the same plot for comparison. 

The results for the basic (unreinforced) panel, 
shown in Figures 21 to 26 of Appendix A, demonstrate that 
the curvatures of the panel, which are of the order ii) ee 
do not have any effect on surface strains; both sides of the 
panel have identical values. In addition, all six figures 
show that the behavior is generally smooth and continuous, 
except near the very edges of the hole for Ee, vs Z/R and 
E.2 VS Z/R. All of the figures also show that the strains 
have virtually achieved their farfield values by 2-4 radii 
from the edge of the hole, which again is similar to isotropic 
plate behavior. 

The computed results for the other configurations 
(shown in Appendices A-G) have several aspects in common 
which are important to understand before discussing those 
aspects which differentiate them. First, because of the 
asymmetry of the reinforcements, the values of Y, distance 
From the reference plane, for the front and back surfaces 
are different. The value of y for the back of the panel 
1s constant at 0.056 inches (for all configurations) and 
the value of y for the front of the panel varies from .112 


to 0.56 inches, depending on the configuration and distance 


from the center of the hole. The large separation of the two 
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surfaces coupled with the large curvatures in the areas of 
greatest reinforcement thickness create a large difference in 
the strain results for the two surfaces. These curvatures 
are created by the reinforcement being pulled toward the 
midplane of the panel as the load is applied, a phenomenon 
well documented to isotropic and orthotropic materials [Ref. 
3]. For example, this curvature of the panel causes the 
unreinforced side of the panel (back) to be in tension, and 
the reinforced side to be in compression, so that the surface 
strains of the back of the panel will always be higher than 
the reinforced side. Additionally, the values of E. plotted 
against X/R on the back side for the reinforced configurations 
are typically higher at the edge than for the unreinforced 
panel, and the reinforced side strains are typically two 
thirds of the unreinforced panel values. The strain values 
for both sides decrease with increasing distance from the 
hole and approach the values for the unreinforced panel from 
above and below at distances far from the hole. Nearly all 
of the configurations (except 6 and 6A) show a marked dip in 
ene distribution of Ee, vs X/R in the vicinity of the outer 
radius of the largest reinforcement layer, similar to the 
decrease in stress found just inside a reinforcement in an 
lsotropic plate [Ref. 6] at a reentrant corner. In fact, in 
nearly every case, the lower portion of the curve is at 

Z2/R = 2.0, a point in between the inner and outermost radii 


meeemost Of the configurations. 
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There are some common features of the results for 
all the configurations which are somewhat unexpected. For 


—€ and e¢ vs Z/R 


le, the curves ¢€ X/R and € 
example, x vs Vf aE _ 


for nearly all the configurations show consistency in general 
character and also in the numerical values close to the edge 
of the hole, in spite of the wide variation of reinforcement 
configuration geometry. The differences in the values plotted 
against Z/R are so slight that they could be ignored to a good 
approximation. Another example is the rather close agreement 
among configurations of the numerical value of the strain 
at the point closest to the hole on the unreinforced (back) 
surface of -the panel (in that location). Every configuration 
has a maximum strain of between 3100 and 330 u strain, values 
in excess of the corresponding value for the unreinforced 
panel. 

Perhaps the most important common feature of all the 
results is the fact that the B subconfigurations (those with 
2 eg layer/layers closest to the basic panel) consistently 
produce significantly lower values of Ee. and Ee, on the 
ieent Surface near the edge of the hole than that for the 
corresponding A subconfigurations. The amount of reduction 
varies among the configurations, with the greatest reductions 
produced by the configurations with larger thicknesses of the 
0° portion of the reinforcement. The DIAL output for the A 
and B subconfigurations consistently shows differences in 


corresponding values of strains and curvatures. The difference 
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in curvature varies with the distance of the elements away 
from the. hole edge, with the greatest differences (at the 
hole edge) being up to 20% for values k, and up to 50% 

for Sa and Koo: In addition, it often happens that corres- 
ponding values of kK. OG ies for the two subconfigurations 
have differing signs which can significantly alter the 
corresponding surface strain values. The midplane strains 
for corresponding locations of the subconfigurations are 
consistently lower for the B configuration than those confi- 
guration A for areas within the reinforcement, and this 
effect often extends out beyond the reinforced region. The 
values of Ee, seem to be most affected, with the variation 
being as much as 20-25% depending upon the configuration. 

In contrast to the common aspects of the results, 
each configuration has some individual results which warrant 
detailed examination. The DIAL output for the first configu- 
ration, shown in Figures 33 to 38 of Appendix B, show several 
interesting features. First, in Figure 33, the difference 
between front and back surface strains (e,) 1s quite large, 
with the corresponding values near the hole edge differing 
by almost 1800 Ho: The maximum front surface strain is about 
1400-1700 Woe depending on subconfiguration, and the minimum 
strain is about 5-700 Ho, OF only 20% of the maximum value of 
the back side. There is little difference in the strains for 
the back surface between the two subconfigurations, showing 
that changes in orientation of the reinforcement material has 


little effect on strains of the unreinforced side. Figure 34, 
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showing Ee, vs X/R, also shows significant difference between 
See subconiigurations and the two surfaces of the panel. 
Figure 35 shows that there is less difference between the 
configurations for Eye! but there is still a large difference 
between the front and back surfaces, on the order of 200 He 
(30% of the maximum value). When Ee, 1s plotted vs Z/R 
(Figure 36), there is a large difference in Ee, between the 
two surfaces and the strain is rapidly changing in the vicinity 
of the hole; however, there is little difference between the 
subconfigurations. Figure 37 and 38 show similar results for 
S and Ey, VS X/R, again showing little variation between 
the subconfigurations. 

The second configuration with a reinforcement, con- 
Sisting of one layer of 1.25 inch radius and one layer of 1.0 
in radius, exhibits behavior which is essentially similar to 
the first configuration, although it differs in degree. In 
Appendix C, Figure 5l plots Ee, vs X/R and shows surface 
strains at the edge of the hole differing by about 1000 u,.. 
The maximum strain on the front side is about 21-2300 He 
and the minimum is approximately 750 Ho: There is far less 
variation between the subconfigurations than for configuration 
ile 

The third configuration results found in Appendix D 
have some slight differences from the two configurations 
above. First, in Figure 61, the maximum front surface strain 


ce) is slightly lower, being between 1800-2000 Ue with the 


minimum about the same as for #2 at 7-800 Mess The dip in €, 
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is more severe as in the previous cases, and there are 
smaller variations between the subconfigurations. [In 
Figure 62, the maximum difference in Ee, near the hole 
edge is about 600 He: There is only a slight difference 
between the front and back surface strains for the B 
subconfigurations (about 20 WH. near the hole edge) while 
there is a large difference between the surface strains 
for the A subconfiguration. 

The fourth configuration results shown in Appendix 
E are similar to the results for the third configuration. 
Figure 81 shows the maximum front surface value of eé, to 
be about 1800-2000 Us with the minimum around 800 Hes 
The dip in ee is shallower and has more gently sloping 
Sides than in the previous results. Figure 82, showing 
Ee. vs X/R shows the same relative difference between front 
and back as the third configuration and a similar difference 
between configurations. Figure 83, showing Ey, vs X/R 
shows the characteristic large difference between surfaces 
for the A and B subconfigurations. 

The fifth configuration results are found in Appendix 
F show some similarity to the first configuration. Figure 
95, showing Ee, vs X/R, shows the maximum front surface strain 
to be about 14-1800 is with the minimum between 500-700 Hos 
The dip in E. is not as deep as the first configuration, 
and it is slightly offset to X/R = 1.5 due to the smaller 
SuW@eer radius. In Figure 96, the variation of ¢€ is seen 


x 


to be about the same as in the other configurations. 
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The sixth configuration, reinforced by only one 
layer, has some major differences from the other configura- 
tions. In Appendix G, Figure 109 shows the maximum front 
surface value of E., to be in the range of 24-2600 Hos with 
essentially no dip to values much below the farfield results. 
This plot also shows the smallest variation between A and B 
of all the configurations, although B still has lower values 
than A. Figure 110 shows that the difference between subcon- 
mrourations for a is not as large as with the other 
configurations. 

To further illustrate the different features of the 
results of the several configurations, the variation of oy 
vs X/R for the front surface of the plain panel and the B 
subconfigurations of the six reinforcements are plotted 
together in Figure 3. This plot clearly shows the range of 
values for the strain at the edge of the holes, the dip at 
X/R = 2.0, and the similarity of values for X/R greater than 
6. To highlight the similarities discussed above for all the 
results, the variation of front surface strains, E.. vs Z/R 
are plotted in Figure 4. This plot shows vividly the 
remarkable consistency between configurations and seems to 
indicate that there is little reduction in strain for any of 
the configurations modeled, despite the variations in height, 
diameter or ply orientation. 

2. Contour Plots 
The output of the analysis computed by DIAL also 


N N 


included contour plots of the Stress Resultants, Noe ee 
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and the out of plane deflections (y direction). The plots 
for the basic (unreinforced) panel shown in Figures 27 to 30 
of Appendix A show that the panel behaves in an isotropic 
manner. The contours all exhibit very localized behavior in 
the vicinity of the hole and comparison of the plots of Ny. 


No and Ny show only one nonzero stress resultant, Nos exists 


Z, 
at an appreciable distance from the hole. At the edge of 
the hole 90° from the applied load (Figure 27), the value of 
N, approaches 3000 lb/in, or three times the applied stress. 
At the edge of the hole at the load axis (Figure 28), the 
value of Ny approaches -1000 1lb/in. Both of the values are 
identical to the values predicted by isotropic theory [Ref. 
7]. Figure 30, the contour plot of out of plane deflections, 
shows that the panel experiences very slight deflections 
along its entire width in a smooth gradual manner. There is 
also a roughly triangular region of deflection which occurs 
above the hole due to the compressive strains in this area. 
Gontour plots of Nye Nye es and out of plane deflec- 
tions are also included in the other Appendices for most 
subconfigurations. Like the strain results discussed above, 
these plots also contain some similarities that are of 
interest. First, as expected, the contours of Ne roughly 
follow the circular shape of the reinforcements near the 
edge of the hole because of the changing thickness in this 
area. The spacing between the contours in the reinforced 


region varies by configuration, being farthest apart where 
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the reinforcements are thickest. Second, the contour plots 


of N, and N, show great similarity among all the subconfigu- 


Zz 


rations, just as the results for Ey and Eye 

The contour plots of out of plane deflections (for 
those cases where plots are available for both subconfigura- 
tions) indicate one important difference between the subcon- 
figurations. With the exception of Configuration 5, all 
other plots show that the B subconfigurations experience 
greater out of plane deflections and that these deflections 
have slightly higher gradients in the vicinity of the hole. 
In addition, the deflections for the B subconfigurations 
seem to extend across the entire panel, while the A subcon- 
figurations show contours of roughly elliptical shape centered 
about the hole. Since all other factors between the subcon- 
figurations are identical, this difference in the nature of 
the deflections must be attributed to the orientation of the 
fibers of the innermost ring. 

3. Computed Strain Concentration Factors 

To reiterate, the results discussed above and shown 
in Appendices are for values at the surfaces of the panels 
above the chosen integration points of the finite elements. 
Because of the established behavior of plates containing 
holes, the quantities of most interest are the values of 
Strain at the hole edge at two positions: 90° and 0° to the 
axis of the uniaxial load. There are empirical and theoreti- 
cal results for the strains at these locations in the litera- 


ture, so it is useful to compare the results of these finite 
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element analyses to those found in the literature. Because 
the empirical results do not specify the precise position 
at which the strain or stress concentration applies, this 
investigation assumed that the point of comparison is the 
midplane of the basic panel. To this end, the output results 
of the finite element program for the elements of the basic 
panel were used as a basis for determining the strains at the 
hole edge. To accomplish this comparison, some means of 
extending the results from the closest point computed 
(X/R, Z/R = 1.05) must be used. The technique that produced 
the most useful results was a simple power curve fit of the 
three values of midplane strain for the panel elements closest 
Memthe hole edge (X/R, Z/R of 1.05, 1.302 and 1.552). The 
curve fit in all cases was excellent when only three points 
were used. When additional points were tried, the accuracy 
of the fit quickly deteriorated because of the generally wavy 
character of the curve near X/R = 2.0, the vicinity of the 
reinforcement edge. Values extrapolated to the hole edge 
using the power curve fit were then used as the numerators 
for the strain concentration factors, Koo and Ko - 

Because the different reinforcement configurations 
Caused the strain distributions to be different, a common 
reference strain, that for the panel without a hole, was 
used as the denominator in calculating the strain concentra- 
tion factors. This permitted unambiguous comparisons between 
the configurations. A panel without a hole, identical in 


dimensions and physical properties but with fewer finite 
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elements, was analyzed by DIAL and the strain for the mid- 
plane of an element in a position corresponding to the hole 
location was obtained. This value of strain (e,) was then 
appropriately factored to represent the strain of a solid 
panel with the same cross sectional area as the panels 
containing the hole. The resulting net section strain was 
Koo and Koy: 


The strain concentration factors computed in this way are 


used as the denominator of the ratio to compute 


listed in Table II. Because finite element models tend to 
be stiffer than actual [Ref. 8], the strain concentration 
factors produced by this technique were expected to be lower 
than the actual factors by roughly the same order as the 
difference in stiffness. 

The results of Table II can be used to gain an under- 
standing of the effect of reinforcement on the strain response 
of the panel. For example, to show the effect of reinforce- 


ment height, the values of K and Koy are plotted against 


90 
the number of reinforcement layers in Figure 5 and 6. If 
configuration (1B) is not considered due to its much larger 
reinforcement volume, Figure 5 indicates that Kyo has a 
Minimum value in the range of two to three layers, suggesting 
that there is an optimum configuration. However, configura- 
tion 1 does exhibit the lowest overall values of Koo Figure 
6 also shows that the number of layers of reinforcement 
influences the compressive strain at the hole edge (0° to the 


load), but the dependence is less clearly defined. In this 


case, the reinforcements which are both tall and wide seem to 
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TABLE If 


COMPUTED STRAIN CONCENTRATION FACTORS 


CONFIGURATION K 


-.812 


=s 75> 
-./17 


-./785 
=-.762 


-./74 
-./12 


-./83 
2702 


te 796 
-./95 


-./88 
-./96 





have the lowest value of Ko - Figure 5 also shows several 
other items of interest. First, it is apparent that all 
subconfiguration B results are lower than for A, and that 1B 
offers the greatest potential reduction in strain, with the 
reduction being approximately 12% of the plain panel value. 
The close grouping of the values for configurations l and 5 
Suggests that increasing the outer reinforcement radius (all 
other parameters being the same) has only a slight effect on 


the reduction of strain. The results for configurations 3 
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and 4 suggest that the effect of the radii of the layers is 
not as great as the number of layers in the reduction of 
strain. As an additional perspective on the results, the 


value of strain concentration factor representing the 


Koq’ 
maximum value of strain in the panel, is plotted versus the 
percentage of hole volume of the reinforcement. This result, 
shown in Figure 7, serves as an approximate measure of the 
efficiency of the reinforcement. This figure shows that all 
of the configurations near 100% of the hole volume have 
virtually the same strain, and that for significant increases 
in weight (configurations 1 and 3) there is only a small 
decrease in strain. The most significant feature of all of 
the above results is the magnitude of the strain reduction 
caused by the reinforcements. The reduction of strain, based 


on the linear computer analysis for a relatively light load 


condition is disappointingly small. 
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III. ALTERNATE THEORETICAL ANALYSIS 


A. UNREINFORCED PANEL 

Because of the lengthy computation time for even as 
simple as model as described above (up to 1.5 hours of 
clock time on the mainframe used by the author for the 
computation of one subconfiguration) it is useful to 
examine other theoretical analyses to determine if they 
Can produce similar results with less effort. Examination 
of the theoretical solutions for isotropic and orthotropic 
plates containing holes will also put the results obtained 
above in proper perspective. 

Stress concentrations in isotropic plates containing 
holes have been understood for some time. Timoshenko 
[Ref. 9], extending the work of Kirsch, shows the maximum 
stress concentration factor, K, to be 3.0 at a location at 
the edge of the hole 90° from the applied load. The plot of 
K versus distance from the hole shows the rapid decrease of 
stress away from the hole. In addition, K is shown to be 
-1.0 at the edge of the hole at the axis of the load applica- 
tion. Rowlands, et.al. [Ref. 1], show that for a plate 
constructed of orthotropic laminae with the same dimensional 
criteria, the stress concentration factors are dependent on 
the material properties of the laminate. The dependence is 


of the form 
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Using the material properties for the basic laminae and for 
the laminate (see Table V and Table VI), the equations above 
produce the following results 

K = 4.106 


90 


K = -.9845 


0 
Comparing these results to the theoretical isotropic results 
demonstrates the quasi-isotropic nature of the symmetric 


layup. 


B. REINFORCED PANELS 
1. Empirical Isotropic Approaches 

There are few analytical results for asymmetrically 
reinforced isotropic panels, but some experimentally derived 
values of stress concentration factors can be found in the 
literature. Peterson [Ref. 10] shows results for an unsym- 
metric circular reinforcement with B/a (reinforcement 
diameter/hole diameter) and h/t (reinforcement panel thick- 
ness/panel thickness) as parameters. These results, for the 
case where reinforcement volume is 100% of the removed hole 
volume and the ratio a/t equal to .1883, show the value of 


Koo Varying between 2 and 3.0. 
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Since the panel being modeled and tested is quasi- 
isotropic, it is of interest to determine whether or not 
these results can be extended to the panel under investiga- 
tion, even though the panels do not precisely match all of 
Peterson's parameters. If the results obtained by this 
approach are sufficiently close to the results of the finite 
analysis and are verified by testing, using these results as 
a first estimate of the actual behavior of the panel would 
seem to be not only attractive, but justified. 

To adapt the panel reinforcement cross-sections to 
the rectangular form used for the isotropic results, a 
Simplifying assumption must first be made. Since the actual 
reinforcements are layered in steps, the volume of the actual 
reinforcement is adjusted mathematically to an equivalent 
rectangular cross-section with an averaged radius and identi- 


cal volume. With the volume of the reinforcement given by 
= Ca... 50 
V ) gi Seater xs) 


the averaged radius is 


V 


eo gn 
i a) t. 


These equations are used on configurations 1-6 to determine 


the parameters B/a and h/t to determine K from Figure 97 


90 
of reference 7. Using the three net section curves shown 
to extrapolate a Similar curve for a/w = 0.1, values of Kg 


were taken from the graph and listed in Table III. 
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There is another set of results for isotropic 


material which has ability to account for the roughly 
triangular cross sectional shape of the reinforcements 
[Ref. ll]. These results, predictions of gross section 


stress (strain) concentration at the hole edge, are for 


elliptic holes in an isotropic plate for material with a 
Poisson's Ratio of .33, quite close to the value computed 
for the laminate. The results are for symmetric reinforce- 
ments, but can be used for asymmetric reinforcements by 


adjusting the cross sectional area by a factor dependent 
upon its shape. The specified value of the correction 


factor for the roughly triangular shape of most reinforce- 


ments is 1/3. The cross sectional area is used along with 
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the values of major and minor axes of the elliptical hole 
(for this case, a = b = 0.5 in) and the thickness of the 


panel to form the parameter 


A 
(a + b)t 


which is the required known value to find the value of K 
from Figure l of Ref. 10. The cross sectional area A for 


use in this approach is given by 


i ele oan = r,s) 


for the configurations under study. Table IV shows the 
values of the parameter A/(a+b)t calculated for the con- 
figurations and the values of K obtained from the plot of 
Ref. 10 and corrected to net section values by the appro- 
priate factor, i.e., the ratio of gross and net panel cross 
sectional areas. These values are generally lower than 


those produced by the previous isotropic prediction. 


TABLE IV 
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IV. EBXPERIMENTAL METHOD 


A. PANEL SELECTION AND CONSTRUCTION 


The particular layups for the basic panel and reinforce- 
ment configurations analyzed above were dictated in part by 
considerations involved in the experimental portion of this 
investigation. In addition, as previously noted, quasi- 
isotropic panels are the most frequent topics of other similar 
investigations in the literature. The use of such panels in 
this investigation permitted comparison of the results with 
these other investigations. The use of a quasi-isotropic 
layup also has the following advantages: 


1) residual stress from curing are minimized, if not 
eliminated. 


2) the layup minimizes the out of plane deflections 
under loading of the basic panel, which reduces the 
chance of masking deflections produced by the assym- 
metric reinforcements to be tested. 
The particular choice of an eight layer [0/+45/-45/9Q) J. 
layup was the result of two considerations. First, this 
is a universally employed and easily duplicated layup. 
Second, using eight layers produces a laminate of significant 
strength. This high strength permits high test loads, which 
in turn reduce the magnitude of potential experimental errors, 
such as small load fluctuations, in relation to the total 


MOoad. This reduction of the error fraction has obvious 


benefits for the accuracy of the experiment. 
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Because of the limited scope and resources of this inves- 
tigation, the number of panels selected for testing was 
restricted to five configurations. First, to establish a 
baseline for strain behavior and ultimate strength for the 
laminate, a basic panel without a hole was used. This | 
configuration, tested twice for reliability, was also the 
means used to verify that the test apparatus applied an 
evenly distributed load across the panel. The second 
configuration was the basic panel with a hole, also tested 
twice for reliability. This configuration was the baseline 
for strain behavior in the vicinity of the hole for the 
unreinforced configuration. 

The last three configurations, each single tests, were 
three of the reinforcement configurations analyzed in the 
Computational Analysis section above. Rather than concentrate 
only on those configurations which achieved the lowest strain 
concentration factors, the experimental investigation con- 
Sidered the full spectrum of the results of the analysis by 
testing configurations 1B, 3B and 6B. By examining the 
response range of these three configurations and comparing 
the experimental results to the corresponding finite element 
predictions, more information could be gained about the 
analysis and the potential effects of the reinforcement 
parameters on the computed results. These three configura- 
tions were chosen because they represented the widest spread 


of the number of reinforcement layers, outer radius of the 
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reinforcement and volume of the reinforcement. All of the 
test panels were tested to failure to reveal potential non- 
linear effects or degradation of performance caused by 
bending of the reinforcements at loads above the analysis 
load. 

The physical dimensions of the panel (shown in Figure 8) 
were the result of two considerations: limitations placed 
on the panel size by the test apparatus, and requirements 
dictated by the choice of hole diameter. Because of the 
well known variation of strain concentration factor with 
hole size [Ref. 12], a diameter of 1 inch was used to avoid 
these effects. The chosen diameter is also representative 
of the potential applications mentioned previously. With 
this chosen diameter and the first estimates of reinforce- 
ment diameter, the 10-inch width of the panel provided more 
than the 5 to 1 panel width to hole diameter ratio suggested 
by Timoshenko [Ref. 8]. The twenty-inch untabbed length of 
the panel was intended to provide sufficient distance for 
uniform load distribution. The width of the tabbing provided 
sufficient distance for load diffusion into the panel. 

The material properties of the HMF 330/34 graphite epoxy 
cloth used for the preliminary design calculations and the 
finite element analysis of the panels are listed in Table V. 
In addition to these values of laminae strength, carpet plot 
data for HMF 330/34 laminate properties are also available. 


Table VI shows the values for the basic quasi-isotropic 
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laminate. The ultimate strength predicted by similar 
carpet plot data for uniaxial loading 1S approximately 
60,000 psi, which for the basic panel corresponds to a 
tensile load of approximately 67,000 pounds. All data 
used in the tables below were furnished by Lockheed 


Missiles and Space Co., which manufactured the panels 


tested. 


TABLE V 
MATERIAL PROPERTIES OF HMF 330/34 


ROOM TEMPERATURE DATA NOT COMPENSATED FOR MOISTURE DEGRADATION 


TENS ION 

ote 11.1 msi E7? 10.4 msi 
COMPRESSION 

a 9.5 msi nA? 8.9 msl 
SHEAR 

a 1.Q msi 

interlaminar shear 150 ksi 


ADDITIONAL PROPERTIES 


att 1e6 uin/in/°OF 
22 20 uin/in?F 
density .056 1lb/in? 


Mi, 





TABLE VI 
QUASI-ISOTROPIC LAMINATE PROPERTIES 


For 503 +45° Plies 


oo 8.15 x 10° psi 12 ae 
aoe 7 ogee 10° psl 
ene 2 ek 10° psi 


B. TEST APPARATUS 

A Riehle PS 30Q Test Machine was used to apply the 
uniaxial tensile loading to the experimental panels. The 
concentrated applied load was diffused through the panel 
by a whiffle tree, shown in Figure 9. Because of the large 
load needed to test the panels to failure, it was necessary 
to apply the load through the four bolts at the upper edge 
of the plates and diffuse the load through the plates by 
the eight lower bolts which pass through the tabbing of the 
panel. The eight lower bolts were torqued to a moderate 
value (90 in-lbs) and the faces of the plate which mate to 
the fiberglass tabbing were checkered to insure that the load 
diffused completely into the tabbing and the panel. The 
whiffle tree was designed to withstand a load 1.5 times 
greater than the anticipated failure load of approximately 
67,000 pounds for the basic panel without hole or reinforce- 
ment. All parts of the whiffle tree were constructed of 4130 


steel and all bolts were close tolerance high strength bolts. 
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Strain gauges were used to record the response of the 
panel to the loading. Other methods such as photoelastic 
coating techniques were not used because of the risk that 
the slight reinforcement from those materials could change 
the actual behavior of the asymmetric reinforcement from 
that which would occur in the absence of the coating. In 
addition, because of the nature of the panels to be tested, 
only the unreinforced side of the panel could have been 
coated. This would then have had to be considered in the 
analysis of the panels. This additional effort in the 
analysis was thought to be an unnecessary complication. 

Because of the rapidly changing strain states near the 
edge of the hole, small single element gauges (120 ohm) were 
used to enable close grouping. A total of twenty-one gauges, 
including 3 gauges in one rectangular rosette, were arranged 
as shown in Figure 10. Because of the symmetrical nature of 
the loading of the panel, all gauges on the same side of the 
hole were oriented in the same direction and different sides 
of the hole had different gauge orientations. In addition, 
gauges were mounted tangentially to the hole edge at different 
orientations to the load axis to investigate the behavior of 
Strains at those locations. Three gauges were placed on the 
back of the panel at the most likely locations of out of 
plane deflections, at orientations of G2 45° and 907 to 


the load axis (gauges 19, 20, 21). 
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NOTE: Gages 19, 20, 21 directed behind gauges 
>, 15, 9 on unreinforced Surface. 


Figure 10. Strain Gauge Placement On All Panels. 
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The strain gauge measurements were made using a Vishay 
Measurements Group System 4000. This system consists of a 
controller, 5 Measurements Group Model 4270 Strain Gauge 
Scanners mounted in parallel (100 channel capability), a 
Hewlett Packard 9825T desktop computer and software to 
operate the system. The software is capable of assigning 
a channel for automatic temperature compensation, configuring 
any channel for differing types of strain gauges or sensors, 
incorporating calibration data, gauge factor, transverse 
sensitivity coefficient, and material property definition. 
The system records and reduces all zero readings and all 
subsequent readings which can be manually or automatically 
initiated. For this experiment, only manually initiated load 


readings were used as discussed in the section below. 


fee LEST PROCEDURE 
To obtain the experimental data for each of the panels 
tested, the following steps were used: 
1. With the panel attached to the whiffle tree and 
hanging in the Riehle machine, strain gauges connected 
to the scanners and system 4000 up and running, zero 
readings were taken. 
2. The whiffle tree was then clamped into the lower 
set of jaws and the panel was subjected to small tensile 
loads (1000 lbs) for several cycles until the strain 


readings for two cycles were virtually the same. 
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3. The lower jaws were released and a new set of zero 
readings were taken and recorded. 

4. The lower jaws of the Reihle were reclamped onto 
the test item. 

5. The load was raised to and held at predetermined 
levels and strain readings were recorded for later 
reduction. 

6. At several points, strain readings were printed out. 
Based on these results, the subsequent load increments 
could be changed to obtain more readings before the 
anticipated failure. 

7. As the loading progressed, significant events such 
as fiber cracking and noticeable deflections were 
recorded in writing and photographically. 

8. Testing proceeded until total failure of the panel 
was achieved. 

9. All collected data were reduced and printed; 


additional photographs of the failed panel were obtained. 
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V. EXPERIMENTAL RESULTS 


A. UNREINFORCED PANEL 

The experimental results for the two panels without 
holes established the baseline behavior of the laminate 
without reinforcement. These results confirmed that the 
Strain response of the panel is linear up to the point of 
failure, with an average failure load of 65,000 lbs. The 
strain (e.) for a 10,000 lb load was recorded and corres- 
Bemas to 1515 Ue for the net section of the panels with 
holes. This is approximately 18% higher than the value for 
panel strain based on the DIAL output. 

The two basic panels with holes were tested to failure 
and produced consistent results. Their ultimate tensile 
loads varied slightly, with panel l failing at 35,000 lbs 
(34,700 psi based on net section) and panel 2 failing at 
37,000 lbs (36,700 psi). Failure was preceded by audible 
Fiber cracking beginning at approximately 18,000 lbs with 
increasing intensity and frequency up to total failure. 
Both panels failed very rapidly, with no visible deformation 
Or damage prior to failure. The edge of the hole at the 
fracture showed signs of delamination just inside both 
surfaces of both panels. The panels failed at the edge of 
the hole 90° from the load through the horizontal axis of 
symmetry, with the fracture continuing in a nearly straight 


line to the side edge. The response of the strain gauges for 
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both panels was linear up to the point of failure, and 
corresponding gauge readings for the two panels showed 
Mistgnificant variation. Figures ll, 12 and 13 show the 
experimental measurements and computational predictions 

EOr Ee. and Ee, vs X/R and Z2/R for a 10,000 lb tensile —— 
load. The values shown are the actual strains plotted at 
the positions corresponding to the location of the middle 

of the gauge. In general, the experimental values show 
excellent agreement with the analysis, particularly close 

to the edge of the hole where the geometric size of the 
elements was smallest. To underscore the linearity of the 
strain response and the absence of significant warping of 
the panels, Figure 14 shows the strain readings for one pair 
of corresponding gauges (5 and 19) on the front and back 
surfaces of the panel 90° from the load axis. Similar 


gauge pairs at locations 45° and 0° from the load axis 


produced similar results (Appendix A, Figures 3l and 32). 


B. REINFORCED PANELS 

The reinforced panels which were tested (configurations 
1B, 3B and 6B) had several aspects in common with the basic 
panel results. During the tests, none of the reinforced 
panels showed observable deformation or visible damage prior 
to failure. The reinforced panels also had the characteristic 
fiber cracking noises. The onset of audible fiber cracking 
appeared to be delayed slightly, beginning at loads approxi- 


mately 21,000 to 23,000 lbs, depending on the configuration. 


65 





(MICROSTRAIN) 


Sy Iaith 


“OOO. 


Ao 
e000 . 
reOG. 
ey. . 
oe 
Baiguane. 1% 


Comparison of Computed and Measured 
vs X/R Basic Panel. 


Strains 


S 
Z 


66 


COMPUTED 


PANEL #1 


PANEL #2 


X7R 


ee 


A 
a2 





-—~ 
poe 





(MICROSTRAIN) 


Sy se Valuilan 


me Si). 


i106 


eo . 


OU C. 


me 0s 


Figure 12. 





COMPUTED —_— 

PANEL #1 (ES 

PANEL #2 4 
2 Y G. 8 


Comparison of Computed & Measured 
Strains Ee, vs X/R - Basic Panel. 


67 





eg In 


20. 


iS. COMPUTED —— 
PANEL #1 © | 


PANEL #2 LA 


ee 
© 
os 
O's 
oe 
2 eee eee 
Q : | ! 
~1200. -i0Q0Q0. -3800. -600. -UQOQ. -2au. ay 


STRAIN (MICROSTRAIN) 


Figure 13. Comparison of Computed and Measured 
Strains E, vs Z/R - Basic Panel. 


68 





Wiis Wisi Clea) (ele ls) 


Or eavel is 


Ldo0vd. 


@ 
Q 
® 
A 
© 
A 
by © 
8O0G0. o4 
‘oho 
A 
© 
A 
© 
600d. = 
© 
A 
uggd. 
g © FRONT (GAUGE 5) 


ZA. BACK (GAUGE 19) 


2 OOO 


Q. : 
ie 8 6.184)" 21000) sCue UF 400900 


LGAD (LBS) 


' O 
Figure 14. Comparison of Front and Back Strain 90 
From Load Axis - Basic Panel 
Gauges 5 and 19. 


69 





All of the reinforced panels failed virtually instantaneously 
and the fracture appearance in all cases failed identical to 
that of the basic panel. A typical fracture is shown in 
Figure 15. In addition, all reinforced panels showed evidence 
of delamination in the region of the fracture, both at the 
edge of the hole and at the side edges of the panels. The 
delamination extended roughly 10° either side of the 90° 
position and was most severe in the two layers closest to the 
front and back surfaces. The comparisons of measured and 
computed strains for a 10,000 lb tensile load (Appendix B, 
Figure 48 to 50; Appendix D, Figures 75 to 77; Appendix G, 
Figures 123 to 125) show the same agreement seen in the basic 
panel results, with the closest agreement occurring near the 
hole edge. 

There are some results for the reinforced panels which 
have both similarities to and differences from the basic panel 
results. One example of this 1s the behavior of strains 


Oo Oo 


measured tangentially at the hole edge One S04. 45, 60° 


and 
90° to the axis of the load. The tangential strains measured 
for the basic panel and reinforced configurations for a load 
of 30,000 lbs are shown in Figure 16. The basic panel strains 
show a smooth transition from compressive to tensile strains, 
very similar to results in section 2.5.1 of [Ref. 13]. The 
unreinforced panel has the highest value of tensile strength 
at the edge of the hole 90° from the load. This plot also 


Shows that the strains for all configurations are very similar 
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Figure 15. Typical Panel Fracture. 
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Figure 16. Measured Tangential Strains at 
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at the top of the hole (0° to 30°) but they begin to separate 
in the region of increasing tensile strain. In addition, the 
reinforced panel strains are reasonably close to the basic 
Panel strains at the top of the hole, and then increase more 
slowly (depending on the configuration) than the basic 
strains, with the widest separation between all the results 
Securring at the 90° location. 

The added thickness of the reinforced panels permitted 
the inclusion of a gauge mounted in the hole at the edge 90° 
from the load axis, to measure e, at the exact hole edge. 
The strain measured by this gauge was included in the compari- 
sons to the computed strains (see plots of E, vs Z/R in the 
Appendices) and in all cases appears to be close to the 
average of the two surface strains (extended to Z/R=l). To 
study the behavior of strain at the hole edge, the measured 
strains for all three configurations are shown plotted 
against the applied load in Figure 17. The plot indicates 
that there is a slight difference in stiffness between the 
configurations and that all the panels behave nearly linearly. 
At loads above 27,000 lbs, Configurations 1B and 3B show a 
slight increase in slope of the data points, indicating that 
some stress relief was occurring around the hole. This stress 
relief, described in [Ref. 14], is caused by microfailures of 
fibers, usually in regions of maximum strain. In addition, 
both configurations show values of strain in excess of 10,000 
Me (or 1% strain), which is the generally accepted limit of 
strain for this material. 
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There are some results for the individual reinforced 
panel configurations which are unusual and important. The 
panel reinforced by Configuration 1B failed at a load of 
39,500 lbs, an increase of 10% over the basic panel average. 
The surface strain recorded for the back surface (maximum) 
before failure was 46400 Hoe Significantly below 1%, while 
the strain measured in the hole at the edge was greater than 
1%. This could signify that the strain is changing very 
rapidly away from the hole, more so than for the other con- 
figurations. Figures 45 to 47 (Appendix B) show the strains 


O 


on the front and back surfaces at the hole edge 90°, 45” and 


a- 


to the load axis. These plots show major variation with 
the corresponding basic panel plots, and show that the 
region most affected is the region 45° to the load. This 
plot has particularly large differences in values of points 
and the apparent slopes of the curves, compared with the 
other configurations. 

The panel reinforced by Configuration 3B failed at 38,000 
lbs, with a failure strain of approximately 10,000 H (1%) on 
the back surface at the hole edge and over 13,000 He recorded 
by the gauge in the hole. Figure 78 (Appendix D) shows a 
smaller variation in strain between the surfaces 90° to the 
load than Configuration 1B. Figure 79 (Appendix D) for the 
strains 45° to the load again shows a large difference in the 


values and slopes between the surfaces. In addition, there is 


a distinguishable nonlinear region close to the failure load. 
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The. panel reinforced by Configuration 6B also failed at 
38,000 lbs with the back surface strain measured in excess 
of 10,000 H. (the gauge in the hole failed well below the 
failure load). The front and back strain comparisons 
(Figures 126 to 128, Appendix G) show very little difference 
from the corresponding plots for the basic panel. 

To summarize the experimental findings, the relevant 
results for all panels with holes are collected below in 
maple VII. The last column contains the ratio (P/P)) of 
the panel failure load to the failure load for the same 


Sized panel without a hole. 


TABLE VII 


EXPERIMENTAL FINDINGS 











SURFACE 
FAILURE 
STRAIN 







EDGE 
FAILURE 
STRAIN 





FAILURE 
LOAD 


(LB) 
















(gauge 19) 


G 


(gauge 22) 





PANEL 








i 








BASIC #1 10,000 
BASIC #2 10,0007 
1B 10,000 
3B 12,000 
6B 13,000 






NOTE 4. Estimates based on extrapolation of experimental 


data. 
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C. EXPERIMENTAL DETERMINATION OF Koo 
For comparison to the experimental results, the value of 
merain concentration factor, Kog was computed from the experi- 
mental results. The value of strain for the laminate without 
a hole was obtained from the plain laminate tests described 
above, with the value cited on page 64 for strain (e,) used 
as the constant reference strain. Because the basic (unrein- 
forced) panels had no gauges mounted in the hole to measure 
a at the edge, a curve fit analogous to the method used in 
the Computational Analysis section was employed. Using strain 
readings for the first three gauges OT ES along the X axis 
(5-8), at a load of 10,000 lbs, the strain at X/R = 1 was 
computed and used to determine Kag for the basic panel. For 
the reinforced panels, the value of strain recorded by the 
gauge in the hole at the same load was taken to be the strain 


at X/R = 1 for the midplane of the panel. The strain concen- 


tration factors computed with these results are listed in 


fable VIII. 
TABLE VIII 

EXPERIMENTALLY COMPUTED VALUES OF Koq 
Configuration Koo 

Basic 2.42 

1B 1296 

3B 2ac8 

6B 2278 
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Vil OtoetooLONFOF RESULTS 


A. COMPUTATIONAL RESULTS 

The computer analysis for the basic panel and reinforced 
Bemtigurations indicated that the strains and curvatures near 
the hole edge had characteristics which were dependent on the 
configuration geometry and the layer orientation. The output 
from DIAL showed that the values of panel element midplane 
Strain at the Gaussian guadrature point closest to the hole 
edge varied from 2987 He (basic panel) to 2631 vee (panel 1B) 
for - and -1032 He (panel 1A) to -725 He (panel 1B) for Ey: 
The curvature kK, which contributes to Ee, at the surfaces of 
the panel was in all cases for the B subconfigurations 
greater near the hole edge, and decreased uniformly with 
distance from the hole. For the A configurations the value 
of > was typically 15-20% less than that for the B configura- 
tions at the hole edge. While the A configurations showed a 
roughly exponential decrease with distance from the hole, the 
B configurations had a local minimum of k, around eRe P= 13s 
The other curvatures, Ky and Keo were larger in all cases for 
the A subconfigurations than for B, but often included changes 
in sign (sense of the curvature) with increasing distance from 
the hole and for corresponding locations of the subconfigura- 


tions. 
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The increased curvature (k,) fOmetie B subconftigqurations 
would seem to dictate higher maximum surface strains for those 
cases. However, the higher surface strains which would be 
expected do not develop because the midplane strains for the 
B configuration are less than those for the A configurations. 
The lower strain computed for the B subconfigurations have a 
greater number of 0° fibers (the load direction) causing the 
panel to be stiffer with respect to the 0° unidirectional 
tension load applied. This result is consistent with the 
findings of [Ref. 2], which shows that the hole gets ‘softer' 
with additional 45° plies. Therefore, even though the panel 
deflections and curvatures are greater for the B subconfigura- 
tions, the value of Koq 1s based on the strain of the midplane 
of the panel element. 

Because the critical factor during panel testing appeared 
to be the maximum strain (e,) occurring on the back surface, 
the method of calculating the strain concentration factor was 
reconsidered. A second series of calculations were carried 
out which took the DIAL output and computed the strains at the 
midplane of the outer laminae of the back surface (.049 inch 
from the panel centerline) again using Equation 1. Using the 
Same procedure as before, an identical curve fit was used to 
extend the computed results to the edge of the hole for 
comparison to the experimental results. The results for Kgq 
are shown in Table IX. From these results, panel 1B would 


seems to have the highest strain at the back laminae and 
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presumably would fail first, followed by panels 3B and 6B. 
As described in the experimental findings, this is the 


reverse of what actually occurred. 


TABLE IX 


STRAIN CONCENTRATION FACTOR BASED ON MAXIMUM STRAIN 


CONFIGURATION 
1A 
1B 


2A 
2B 


3A 
SB 


4A 
5B 


5A 
5B 


6A 
6B 





To briefy explore the effect of reinforcement shape on 
the results, a computational analysis was completed for a 
one layer square reinforcement to be compared to Configura- 
tion 4. The square patch was sized to have a volume equal 
to that of the hole and had 2.625 inch sides. The patch was 
oriented so that the distance along the X and Z axes was 
1.3125 inches and the sides of the square were perpendicular 


to the axes. Computations identical to those performed in 
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the Computational Analysis section were performed. The 
results indicated that the strains, curvatures and deflections 
Peoetne Square patch and circular patch were virtually identi- 
cal. To verify these findings, a panel with the square patch 
reinforcement was tested. The failure load, strains at 10,000 
lbs and the failure appearance were nearly identical to the 
Circular patch results. These limited results seem to indi- 
cate that the planform shape of the reinforcement has limited 
effect on the strain response of the panel. 

The effect of symmetry was also briefly explored by 
modelling two symmetric circular reinforcements similar to 
configuration 6 for analysis. The first model was identical 
in the material and outer ring radius to configuration 6 with 
a one layer reinforcement 1.5 inches in radius added to both 
sides of the panel, with a reinforcement volume of 200% that 
of the hole. The other model used one layer reinforcements 
on both sides with a radius of 1.125 inches to create a 
reinforcement with a volume equal to that of the hole. The 
same method of analysis was used and the same output generated 
by DIAL. As for the unreinforced panel, the curvatures were 


reduced to the order of 19719 


, With virtually no out of plane 
deflections. The value of Ee, at the edge of the hole (using 
the same curve fit) for the large symmetric reinforcement 

was computed to be 2691 Wor which gave this reinforcement a 
strain concentration factor of 2.09. This represents an 


improvement of 15% over the basic panel and 25% over the 


results for configuration 6B. The smaller reinforcement had 
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Similar results, with a value of Ee, at the hole edge of 
Bo 2 He and Kog of 2.15, representing improvements of 113% 
and 20% over the basic panel and Configuration 6B, 
respectively. 

The results above concerning the variation of reinforce- 
ment volume and planform shape indicate that there may be 
other asymmetric reinforcement configurations which have 
Superior performance to the configurations covered by this 
investigation. The results for Configuration 1B and the 
large symmetric reinforcement clearly indicate that large 
reinforcement volumes do produce greater maximum strain 
reductions. The fact that the square planform reinforcement 
did not have a lower strain concentration factor than the 
comparable round reinforcement implies at the very least that 
reinforcement planform shape does not one harmful effects. 
Confirmation of this planform shape independence could be 
quite useful in applications, permitting the most easily 
producible shapes to be used. However, because only one 
non-circular planform was modelled, it is not possible to 
assert with certainty that the shape of the planform is 


unimportant. 


B. EXPERIMENTAL RESULTS 
The experimental results appear to agree with most of the 


computational predictions regarding strain (e,) near the hole 


edge 90° from the load axis. For the unreinforced panels, 
the difference between the front and back surface strains was 


minimal, as predicted by the computational analysis. The 
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measured strains, Eo averaged 2447 Hos which differs from 

the computed value of approximately 220Q He by nearly the same 
amount (18%) as the difference in net section strains dis- 
cussed above (page 64). Additionally, the value of Se 
from the extrapolated experimental data for the strain at 

the hole edge (3640 Us) is higher than the computationally 
extrapolated value (3176 Ho) by about the same amount. 

For the reinforced panels, experimentally measured surface 
strains (e,) at X/R = 1.25 (location of center of the gauge). 
for a 10,000 lb load are compared to the computed values of 
ec. obtained from the plots of Ee, vs X/R at X/R = 1.25, in 
Table X. In addition, the measured strains from the gauge 
in the hole are compared with the values of Ee, extrapolated 
to the hole edge for the panel midplane in Table XI. From 
these results, Configurations 3B and 6B are in excellent 
agreement with the computed predictions for front and back 
strains, while 1B shows the experimental readings to be 25-30% 
higher than predicted. The difference in readings for 1B may 
be explained by the fact that the finite element model is 
stiffer than the actual panel (previously discussed) and that 
this particular model has three sets of elements layered on 
top of each other near the hole, causing that region to be 
more affected by the analytical approximation. The midplane 
strain comparisons show 1B and 3B to be in good agreement 
with the computational predictions, the error being less than 


the difference in net section values for Ee. discussed above. 
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The experimental results for the gauge in the hole of panel 


mamee Stgniticantly higher than the computer values. 


TABLE X 


COMPARISON OF EXPERIMENTAL AND COMPUTED MICRO-STRAINS 


CONFIGURATION FRONT BACK MIDPLANE 
SURFACE SURFACE 

1B 

Computed 900 2500 2800 

Measured 1194 1821 2919 
3B 

Computed 15a0 30QQ 2884 

Measured 1538 PA FEA) 3452 
6B 

Computed 1800 250Q 2950 

Measured 1876 2700 4204 


The experimental findings indicated that failure occurred 
at the load for which e, reaches approximately 10,000 Hy. (13 
Strain) somewhere through the thickness of the panel and re- 
inforcement, near the hole edge, and approximately 90° from 
the load axis. For the basic panels, the failure strain at 
the edge of the hole on the surfaces was assumed to be the 
same throughout the panel thickness because of the small cur- 
vatures indicated from the analysis. This assumption could 
not be verified experimentally, but there was no visible 
evidence to the contrary. The results for the reinforced 
panels 3B and 6B indicated that the back surfaces of both 
panels at the edge of the hole also reached the strain re- 


quired for failure, while at the same time, the strains 


84 






recorded hy the gauges in the holes indicated strains far 
higher than required occurred before failure. This appears 
to indicate that either the hole is stretching by some 
mechanism, such as the stress relief caused by microfailures 
[Ref. 13] or that the gauges were not recording accurate 
information, or both. The presence of delamination in the 
failure region suggests the possibility that the hole is 
stretching and the gauges may be affected by the complex 
strain state in that vicinity. The findings for panel 1B 
are even less clear in that the maximum surface strains were 
well below 1% strain, but the gauge recording strain in the 
hole did achieve values in this range. Also, the recorded 
strains for the front and back surfaces (Figure 45, Appendix 
B) show far less variation than the analysis predicted 
(Figure 33, Appendix B), which suggests that some mechanism 
is creating an unusual strain distribution. This complicated 
distribution of strain through the thickness of panel 1B is 
not easily explained, but again the presence of delamination 
in the failure region suggests that interlaminar shear is 
both present and contributory to the strain state recorded. 
In addition to the complicated strain behavior discussed 
above, the results of the front and back surface strain 
behavior confirmed that the panels behave in a nonlinear 
fashion near the hole edge when the panel were subjected to 
loads above that used in the computational analysis (10,000 
lbs). Comparison of the ratios of front and back surface 


strains (e,) S01 from the load axis, show that the ratios 
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decrease when the load increases. This suggests that there 
may be some mechanism which 1s reducing the curvature of the 
panel and therefore the difference in strain readings as the 
load is increased. Because the experimental setup had no 
provision for measuring out of plane deflections or curvature, 
there was no way to determine whether the decrease in front 

to back surface strain ratio was accompanied by reductions in 
curvature. The nonlinear behavior of the panel is even more 
evident when Ee. is plotted vs X/R, as done for panel 1B in 
Figure 18. The most striking feature of this plot is the 
increasing ratio of the highest strain (X/R = 1.25) to the 
lowest (X/R = 2.0). For panel 1B, this ratio increases from 
eof at 10,00Q Ibs to 1.73 just prior to failure, an increase 
of about 10%. Similar results for panel 3B (Appendix D, 
Figure 80) and 6B (Appendix G, Figure 129) show less severe 
increases in this ratio, hut there is clearly nonlinear 
behavior at high loads across the entire panel. Because these 
effects are only slightly above variations that could be 
attributed to experimental error, there is no pressing require- 
ment to alter the computational model to perform nonlinear 
computations. 

Before any firm conclusions can he drawn from the experi- 
mental results, some consideration should be given to the 
accuracy of the experiment. Since there are inevitably small 
Variations in panel construction, gauge bonding, and load 


history, it is possible that results for repeated tests would 
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average to higher or lower values for failure load and 
strain. Therefore, it is difficult to assert with certainty 
that panel 1B had the best results or that 6B had the worst. 
However, Since none of the test results indicate improvements 
of more than 10% of the laminate ultimate strength, in spite 
of the variation of reinforcement configuration, there is 
justification for the assertion that none of the configura- 
tions produce significant strengthening of the panel. 

From the discussion of the experimental results above, 
it 1s apparent that the most important measured results were 
the values of Ee, at various locations of the panel. The 
Variation of co with Z/R was also important, showing the 
compressive strains in the panel parallel to the load axis. 
The variation of Ey with X/R and the results of the rosette 
proved to be of little value. The recorded variation of ES 
with X/R provided an added demonstration that the computa- 
tional models were in close agreement with the experimental 
results. The rosette results verified that the principal 
strains were oriented parallel to the load axis away from 


the hole, a result expected for uniaxial loading. 


C. COMPARISON OF PREDICTED AND MEASURED VALUES OF Ko q 


To evaluate the various predictions of K fOr “SUitaDI II ty 


90 
of use, all of the previous results for the basic (unrein- 


forced) panels and Configurations 1B, 3B and 6B are assembled 


in Table XI. Table XI indicates that the results based on 


Figure 97 of Peterson [Ref. 9] are in question, being 25 to 
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30% high and low. The predictions based on this technique, 
shown with the experimental values in FPigure 19, show no 
real consistency or trend. Because of this, the use of such 
predictions seem suspect, even with the results being con- 
servative, i.e., larger than expected from the experimental 
results. 

On the other hand, the alternate isotropic prediction 
method, based on Ref. 10, gives predictions which are in 
reasonable agreement with the experimental results. The 
results closely follow the correct trend as shown in Figure 
20, with only the reinforced panel being far from the predic- 
tion. Because of the general accuracy, the use of a curve 
such as found in Ref. ll seems to be justifiable in cases 
where estimates of behavior of quasi-isotropic plates are 
desired. 

The predictions based on the computational analysis are 
in better agreement with the experimental results than the 
other predictions. The agreement points to the excellent 
modelling capability of the DIAL code, and to the suitability 
of the techniques used in the analysis, i.e., the method of 
extending the computed results to the edge of the hole. The 
slight variations noted for panels 1B and 6B can be partially 
accounted for by experimental errors and difficulties pre- 
viously mentioned. In view of the results, the DIAL code has 
proven to be an accurate means to predict the true behavior 


of a laminate containing a hole. It seems reasonable to 
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assert that similar results can be expected for nearly any 


isotropic balanced laminate to be modelled. 


TABLE XI 


COMPARISON OF MEASURED AND PREDICTED VALUES OF Koo 


2 Ouia 
PREDICTION 
MEASURED 
COMPUTER 


R.Ae.S 


PETERSON 





90 





Be 


3 


ce 


TT 


cre 


eae |) ce ee] mae | | I | | 


Baroiiime: 19, 


Extrapolated 
forevasw = 0.1 


es 


Oe = 


After Peterson, 
A Figure 97. 


NX EXPERIMENTAL 
_——— After Peterson 


ye 1.u ie ee 
H/T 


Comparison of Isotropic Prediction and 
Measured Values of Kgg> 


oul: 






C3 


x DATA 
a COMPUTER PREDICTIONS 






From uose.5. 
REPORT 65004 


\ 


— 


CORRECTED 
FOR 
NE SECTION 


m8 Der 0.4 0.6 


(2) 
me) 


Peo! T 


Figure 20. Comparison of Isotropic Prediction and 


Measured Values of Kgq° 


92 





VII. CONCLUSIONS 


The following conclusions can be drawn from the results 
of this investigation: 

1. The results for the quasi-isotropic panel modelled 
in the finite element analysis indicate that reinforcements 
of circular shape consisting of two or three additional 
layers of material hold the best potential for the reduction 
of maximum strain at the edge of the hole. The improvement 
is not large, with the maximum strains approximately 10% 
lower than for an unreinforced panel with an identical hole. 

2. From the analysis, the reduction of the maximum 
strain seems to be more dependent on the volume of the rein- 
forcement and the number of layers used (height) than on any 
other parameter regarding the physical shape of the reinforce- 
ment. 

3. The computational results indicate that there is a 
greater reduction in maximum strain for configurations which 
have the reinforcement layers closest to the panel oriented 
parallel to the load axis, even though the out of plane de- 
Flections and curvatures are greater for these configurations 
than for configurations with all reinforcement layers oriented 


a5? 


to the load axis. 
4. The finite element analysis code employed in the 
investigation showed excellent general agreement with the 


experimental results. The observed nonlinear nature of 


= 





strains at some points on the panel under high loads is 
small enough (approximately 10%) that it can be neglected 
for the purposes of this investigation. However, the non- 
linear behavior of Ee, at the point on the panel experiencing 
the maximum strain is likely, to influence the failure 
strength of the panel. 

5. The reductions in maximum strain predicted by the 
analysis were corroborated by the experimental results for 
the panels tested. 

6. The experimental results indicate that asymmetric 
circular reinforcements will only increase the strength of 
the panel by 5-10% of the total laminate strength. The 
slight variation of failure loads for the different configu- 
rations tested indicates that the small range of improvement 
in maximum strain indicated in the analysis is true up to 
failure of the panels, and that configuration geometry 
apparently has little effect on the ultimate strength of the 
panels. 

7. The analytical results for the two symmetric configu- 
rations indicate that far greater reductions in strain are 
possible for even very thick symmetric reinforcements. 

8. The use of the Royal Aeronautical Society plot of 
stress (strain) concentration factors for isotropic materials 
to predict strain concentrations for asymmetric reinforcements 
of a quasi-isotropic laminate seems justified by the experi- 


mental results. 
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VIII. RECOMMENDATIONS 


From the experience gained during the course of this 
investigation the following recommendations are submitted: 
1. Further investigations should be conducted to deter- 
mine whether or not Significant improvements in strength and 
Strain reduction can be made by varying the following param- 
eters of asymmetric reinforcements: 
a. Increasing the outer radius of the reinforcement 
beyond X/R = 3.0 
b. Significantly increasing the volume of the 
reinforcement. 
c. Changing the planform shape of the reinforcement. 
dad. Submerging the reinforcement layers below one or 
two layers of the laminate. 
2. In future experiments, strain gauges should be 
arranged to measure Ee, vs X/R and e, vs Z/R on both sides 
of the panel. The measurement of Ee, vs X/R is not of great 
importance. 
3. In future experiments, Strain gauges should be 
mounted tangentially to the hole at several orientations 
to the load axis on both sides of the panel. Also, gauges 
should be mounted in the hole at locations 0° and 90° to 


the load axis. 
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